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The competition between the staggered flux state, or the d-density wave state, and the d-wave 
pairing state is analyzed in two-dimensional t-J model based on the U(l) slave boson mean-field 
theory. Not only staggered flux of spinon but also staggered flux of holon are considered. In this 
formalism, the hopping order parameter of physical electron is described by the product of hopping 
order parameters of spinon and holon. The staggered flux amplitude of electron is the difference of 
staggered flux amplitude of spinon and that of holon. In 7r-flux phase of spinon, staggered fluxes of 
spinon and holon cancel completely and staggered flux order of electron does not exist. However, 
in staggered flux phase of spinon whose staggered flux amplitude is not tt, fluxes does not cancel 
completely and staggered flux amplitude of electron remains. Thus, the phase transition between 
these two phases, 7r-flux phase and staggered flux phase of spinon, becomes a second order transition 
in physical electron picture. The order parameter which characterizes this transition is staggered 
flux order parameter of electron. A mean-field phase diagram is shown. It is proved analytically 
that there is no coexisistence of staggered flux and d-wave pairing. The temperature dependences 
of Fermi surface and excitation gap at (0, n) are shown. These behaviors are consistent with angle- 
resolved photoemission spectroscopy (ARPES) experiments. 
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I. INTRODUCTION 

Recently the staggered flux state [1,2] has been revived 
as a candidate of pseudo-gap phase of high T c super- 
conductors. The staggered flux state is characterized by 
the staggered orbital-current. Time- reversal-symmetry 
is broken if staggered current is the current of physical 
electron. Weak magnetic signals were caught in recent 
experiments of underdoped YBa 2 Cu306+a; (YBCO) by 
neutron scattering [3,4]. Static alternating magnetic sig- 
nals were also caught by muon spin rotation experiment 
in the vortex cores of underdoped YBCO [5] . The stag- 
gered flux state contains a density wave ordering whose 
symmetry is d x 2_ y 2 [6]. Because of this, it is also called 
d-density wave. Thus, there is the gap whose symmetry 



son methods, the U(l) method and the SU(2) method, 
the physical electron operator is described by a prod- 
uct of auxiliary fermions and bosons. In the U(l) slave 
boson theory, the physical electron operator c; CT is de- 
scribed by a product of an auxiliary spin- 1/2 neutral 
fermion operator f %u called spinon and a auxiliary spin- 
less charged boson operator bi called holon; Ci a = bjf ia 

with a constraint b\bi + f} a fio = 1. Here, the re- 
peated spin index a is summed up over the two spin 
states. On the contrary, in the SU(2) slave boson theory 
[19,20], two auxiliary bosons are introduced. The physi- 
cal electron operator ci a is described by a product of an 
auxiliary isospin-1/2 neutral fermion operator fi a and 
two auxiliary spinless charged boson operators, bn,bi2', 
= = {b\ x f ia + &l 2 W/I/)/V2 with three 



is d x 2_ y 2. Although d x 2_ y 2-g&p exists in both the stag- constraints, i^r^w + h\rhi = 0. Here, hf = {b iu b i2 ), 



gered flux state and the d x 2_ y 2-wave pairing state, the 
structures of the fermion excitation in these states are dif- 
ferent. The staggered flux state should have a segment- 
like Fermi-surface as observed in angle-resolved photoe- 
mission spectroscopy (ARPES) experiments [7-14]. On 
the contrary, d x 2_ y 2-w&ve pairing state has a point-like 
Fermi-surface. The absence of specific heat anomaly may 
be thought to deny the staggered-flux phase. However, 
we remark that there are arguments to justify the absence 
[15,16]. The lack of anomaly is either because of finite 
chemical potential [15], or because that this transition 
can be described by 6- vertex model [16]. 

The high-T c superconductors are doped Mott- 
insulators. Essential physics of these materials comes 
from strong Coulomb repulsion. This physics is described 
by t-J model [17,18]. A method that restrict double oc- 
cupancy is the slave boson method. In both the slave bo- 



ces. 



' f: ,), and t = (ri , T2 , T3) are Pauli matri- 



Wen and Lee proposed based on the SU(2) theory 
[20-22] that the pseudo-gap state is a mixture of stag- 
gered flux state and d x 2_ y 2 superconducting state. These 
orders are dynamical rather than static and that the 
time-reversal symmetry is not broken. They discussed 
this in the region where bosons are condensed, (zn , 2^2) 7^ 
0, based on the 0(4) sigma-model description. It is also 
proposed by Lee and Wen [21] that the static staggered- 
flux order can only exists in the vortex core. Short com- 
ings of this theory are that finite size Fermi-surface is not 
obtained naturally and that the condensation of bosons 
is assumed. The staggered current of physical electron 
and that of spinon are not equivalent in the staggered 
flux [1,2] phase at finite temperature because the stag- 
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gered flux phase of spinon exists above the temperature 
of Bosc-condcnsation of holon [19,20,23]. Above the tem- 
perature of Bose-condensation, holon operators cannot 
be treated as classical numbers (c-number). [24] The cor- 
respondence between the electron current and the spinon 
current is not obvious. 

On the contrary, Chakravarty, Laughlin, Morr, and 
Nayak [25] proposed that the staggered flux order exists 
as a static form in the pseudo-gap state. They proposed 
this at zero temperature based on Ginzburg-Landau the- 
ory. The staggered flux state they proposed is that of 
electron and time-reversal symmetry is broken. They 
also proposed that there exist three ground states. The 
ground state changes as doping increases from pure stag- 
gered flux phase to the coexisting phase of staggered flux 
state and d-wave superconductivity(dSC) state and fi- 
nally to pure dSC phase. A short comings of their the- 
ory is that it is a phenomenology without microscopic 
foundation. 

In this paper, we analyze the staggered-flux state of 
physical electron in two-dimensional t-J model by the 
U(l) slave boson method without assuming the Bosc 
condensation of holons. The relation between the stag- 
gered current of electron and the staggered current of 
spinon is provided. In Sec. II, our formulation will be 
reviewed. Not only an auxiliary field Xij that describes 
spinon-hopping but also an auxiliary field B^ that de- 
scribes holon-hopping will be introduced to decouple the 
hopping Hamiltonian. The advantages of this formula- 
tion will be stated in Sec. II A. For mean field solution, 
not only staggered flux of spinon but also staggered flux 
of holon will be considered. This solution provides the 
relation between the staggered flux state of spinon and 
that of electron. In Sec. Ill, The self-consistent equa- 
tions will be derived. In Sec. IV, the phase diagram is 
shown, where the staggered-flux phase of electron exists. 
In Appendix, it is proved analytically that the staggered 
flux and the d-wave pairing do not coexist. The insta- 
bility of these two states, the staggered flux state and 
the d-wave pairing state, are discussed. In Sec. V, the 
temperature dependences of Fermi surface and excita- 
tion gap at (0, 7r) are shown. We show that the Fermi- 
surface in the staggered-flux state resembles segment-like 
Fermi-surface observed by ARPES experiments. Finally, 
conclusion will be given in Sec. VI. Part of the present 
results will be published elsewhere [26]. 



II. FORMALISM 

A. Slave boson t-J model and auxiliary fields 

The starting Hamiltonian is two-dimensional t-J model 
on a square lattice [18], 

H= -tY J P(clc ja + h.c.)P + jY J Si-S j . (1) 



Here, {i,j) represents sum over the nearest-neighbor 
sites, the repeated spin index a is summed up over the 
two spin states, and P is a projection operator to 
no doubly occupied state, Si — \^ ia (<j) aa > c icr < , where 
<x= (<ti, £72, (73 ) are Pauli matrices. One way to represent 
the projection to no double occupancy is the slave boson 
method. In the U(l) slave boson method, the electron op- 
erator a a is described by a product of an auxiliary spin- 
1/2 neutral fermion operator fi a called spinon and an 
auxiliary spinless charged boson operator bi called holon; 
Cia = b\fiv with a constraint b\bi + f- a fiv = 1. Here, 
the repeated spin index a is summed up over the two 
spin states. In the path integral formalism, the partition 
function Z(0) — Trexp(— 0H) is described by following 
functional integral, Z = J[db][df][d\]ex.p(—S), 

S = f ' driy^ibldrbi + fldrM + H], (2) 
Jo 

II 'E (fUi^bi + c.c.) • •/ E S - ■ 

+i E + f&u, - 1) - Me E f&™ ( 3 ) 

i i 

where is the inverse temperature = 1/T, r is the 
imaginary time, J[db] = J Yli T db\{T)dbi(r) is a complex 
boson integral, J[df] = J l\ iTa . df} a (T)df ia (T) is a com- 
plex Grassmann integral, J[dX] = JY[ir^i( T ) * s the 
Lagrangc-multiplier integral that represents a constraint 

h % + flfi* = 1, Si = yl(cr) aa >f ia >, Me is chemicsl 
potential. We introduce three complex auxiliary fields 
Xij, rjij, and B^ on links to decouple the Hamiltonian. 

<ij> 

+ t{Bi jX * j +B* jX i j } 

^~{Xijfj a fi<7 + Vij.fj(rficr' e aa' } + c - c - 

+Y fel 2 + M 2 }] 

+1 E \{b\bi + flfia - 1) - Me E fUi*- ( 4 ) 

i i 

Then partition function Z{0) is rewritten as; 
Z = J[db][df][d\][dx][dr]][dB}exp(-S'), where S' = 

Ip dT[ZMd T bi + fldrfir) + H] . The following rela- 
tions are obtained by differentiating the integrand in Eq. 
(4) by xtpVij, and 5* [27], 

Xij — (Xij) — (ficrfjv)i 

Vij = (vij) = {fnfji - fnfjx), 

B^ ee (Bij) = (bib,). (5) 

According to these relations, X ij describes the hopping of 
spinon, rjij describes the singlet pairing of spinon which is 
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called resonating- valence bond (RVB), and B^ describes 
the hopping of holon. The quantity Bjj = (b\bj) can 
have a finite value even if there is no Bose-condensation 
of holon(see Sec. IV). 

Present treatment is different from previous one by 
Ubbens and Lee [23]. In Ref. [23], Bij was not intro- 
duced. Then (Xij) can not be interpreted as spinon hop- 
ping {flu f jo), but the correct relation becomes (Xij) — 
(fiafja + §jb\bj) as can be known from Eq. (4) without 
term. In their method, the four-boson term blbjbjbi 
is created by the decoupling, but it is neglected based on 
an argument that the effect of this term is small at low 
doping. The situation is similar in the SU(2) method 
[19,20], if B^ is not introduced. Our treatment has a 
merit that the four-boson term does not appear in the 
decoupling. 



B. Saddle point solution 

We approximate the integral by auxiliary fields Xiji 
rjij, and B^ with their saddle point values Xiji Vij> anQl 
Bij . For the saddle point solution of Xij and B^ , we con- 
sidered not only the staggered flux order of the spinon 4> s 
but also the staggered flux order of the holon <^ s (Fig. 1), 
namely 



X l+SlI = xe l( - 1)>s/4 =x s + i(-iyy s , 
X t+ y, t =xe- l{ ~ iy ^ /4 =x s -i(-l) i y s , 
/>',.,, = BeWM* = x h + i(-iyy h , 



B i+§ti = Be-*(-DVh/4 = Xh _ i(_i)< 



(6) 
(7) 
(8) 
(9) 



Here, x and y are unit vectors in the x and y direction, 
x s = xcos(</> s /4), y s = xsinO s /4), x h = Bcos(0 h /4), 
yh = ^sin(0h/4) [28]. The staggered flux state contains 
a density wave (particle-hole pairing) ordering whose 
symmetry is d x 2_ y 2, which is called "d-density wave" 
[6]. The order parameters y s and j/h correspond to the 
d-density wave order parameter of spinon and holon, re- 
spectively [29]. For the symmetry of spinon pairing, we 



considered d x 2_ y 2 [30,31], namely 



(10) 



N/ <])s /V tya\f \/ (])h/\- ())hN/ 



<- 



The integral by Lagrange multiplier field Xi is also ap- 
proximated with its saddle point value A^ = —iXq. Then 
spinon and holon have chemical potentials, fi — /i — Ao 
and /ih = — Ao, respectively. These chemical potentials 
enforce the global constraints of spinon number and holon 
number, (fjfi) — 1 — 5 and {b\bi) = 5, respectively. 

This method that treats the spinon and the holon as 
separated degrees of freedoms is justified because the low 
temperature phase of the t-J model is in the deconfinc- 
ment phase of spinon and holon even when full gauge 
fluctuations (fluctuations around this saddle point) are 
included. It is shown based on the compact U(l) treat- 
ments of the gauge fields [32,33]. As the gauge field cou- 
ples with finite density fermion, the gauge coupling gets 
weaker. In the deconfinement phase of spinon and holon, 
the gauge field can be treated perturbatively. 

Our formalism is an extension of the study by Ubbens 
and Lee [23]. The term couples not only to 

(3J/8)xij but also to tBij, namely the spinon feels 
the spinon staggered-fiux(0 s ) and the holon staggered- 
flux(0h)- The expectation values of the holon, B and 
4>h, have finite values for the solution of self-consistency 
equations. Two advantages exist in our formalism; 1) 
this is a new saddle point solution whose free energy is 
lower than the previous one, 2) this solution provides the 
relation between the staggered- flux of electron and that 
of spinon. 

With holon order parameters B^, the hopping order 
parameter of electron becomes a product of the hopping 
order parameters of spinon and holon, 



(flfi*){bfa) 



x Be 



±i(0 8 -0h)/4 



(11) 



The electron staggered-flux order parameter <p e and the 
electron d-density wave order parameter y e are given by 



bh i Ve = xBsin(4> e /4). 



(12) 



According to eq.(12), the d-density wave order pa- 
rameter of electron y e can also be written as, y c = 
~*2lv Efc( cosfc z ~ cos/c a )(cj t(T c k+ Q (T ). 



III. SELF CONSISTENT EQUATIONS 

This saddle point is described by the following par- 
tition function Z MF = J[db][df]exp(-S MF ), where 

S MF - Jo" dr^Mdrh + fld T fi„) + H*^] , 



H 



MF 



Hotter ~ NQjl - W.) 

r 3J 



+ 2N (z' + Vl + V 2 ) + M{x s x h + y s y h )] , (13) 



FIG. 1. Staggered flux of spinon <f> s and holon 0h. 



the spinon chemical potential is /i = /i — Ao, and the 
holon chemical potential is /ih = — Ao- The saddle point 
value of iXi is iA = Ao. Here, the matter part of this 
Hamiltonian is given by 
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ttMF 
^matter 



= ^'^A + *W^. (14) 



By minimizing the free energy F, we obtain the self- 
consistency equations, 



where, ^ k stands for a sum over the half Brillouin 
zone. <]/ k and ^{j are vector representations of spinon 
and holon, = (fl p fl +Qp /-k-Qi), *£ f = 

(^k'^iL+Q)- The matrices i? k and are, 



ffk = 



/ ( £ k - m) 
-*Wk 
A k 




-(ek + M) 
-(ck-M) 
-A k -*W k 




(15) 



x tanh(/3S ks /2) 
y s = (2i?/ h + — y, 



2£ k 
3^ 
4 



k,s 



7^- 



-(4 + M h ) 



k,s 



(16) 



where 



tanh(/?g k ,/2) r 



s(-/z) 



(4 + ^ k 2 ) 1/2 



£k 


3J 

= — (2i.T h + — x s )(cos fca, +cosfc y ), 


(17) 


w k 


= (2iy h + ^-y s )(cos fc x - cos k y ), 


(18) 


Ak 


= (cosk x — cos k y ), 


(19) 




e k = — 2ix s (cos k x + cos fcy), 


(20) 




W k = 2ty s (cos k x — cos fc y ). 


(21) 



V = 



3 J 1 v^, 



2 tanh(/3£; ks /2) 
7k- 



k,s 



2£ ks 



(-*) 



2(^2 + ^2)1/2 e /3EL_l' 



2/h 



2 ^E' 

k.s 



2(^2 + ^2)1/2 e ^ 3 _r 



(26) 

(27) 
(28) 

(29) 

(30) 



The spinon couples not only to spinon order parame- 
ters, x s and y a , but also to holon order parameters, x h 
and 7/h- The holon couples only to spinon order param- 
eters, x s and y s . There exist two components, Wk and 
A k , that generate the d x 2_ y 2 wave gap. 

After the Bogoliubov transformation, diagonalized 
Hamiltonian is obtained, 



The chemical potential jj, and /Zh are determined by 



MF 



H, 



where 



E 'EU^L^s-^s) + Ei s allat (22) § = l_ , / /2 _ x 



tanh(/3£ ks /2) 



^ = ^ (syfel + W* A*) 2 + A k , (23) 

^ = s\f¥+WF-f> ( 24 ) 

Both aks and /3 ks are fermionic fields, a ks is bosonic field, 
and E^J is excitation spectrum of spinon (holon) . Index 
s describes the band index, which takes the value +1 or 
— 1. After integrating out spinon field and holon field, 
the free energy has the following form, 



N e 0E ±* - 1 ' 

k,s 



, (31) 
(32) 



F = -2T^'lncosh(/3£; ks /2) 



k.s 



+T^'ln(l-e- 

k,s 

3 J 

+ 2i v[_-( a; 2 + yl + ^2) + 2t( XaXh + y s y h )] 
-N6(n-v h ). (25) 



where 7 k ± = cos k x ± cos 

At half-filling, where spinon chemical potential /i is 
zero and holon is absent (Bij = 0, i.e. Xh = yn = 
0), E ks = ^e k + (3J/4) 2 (y 2 + V 2 )(cosk x - cosfc y ) 2 , and 
both y s and 77 have the same self-consistent equations. 
Thus, staggered flux state and d-wave pairing state are 
degenerate. All states that have the same value of yl + rj 2 
are degenerate. However, which state has the lower en- 
ergy at finite doping is known only after actual minimiza- 
tion of the free energy is done. 
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IV. PHASE DIAGRAM AND STAGGERED 
CURRENT 

A. Phase diagram 




0.02 0.04 0.06 0.08 0.1 0.12 



5 t/J 

FIG. 2. Mean field phase diagram for t/J — 1, where S is 
hole concentration. The staggered-flux order of electron ex- 
ists only in region 2. In this phase, staggered-fluxes of spinon 
and holon do not cancel completely. Thus the staggered flux 
of electron remains. In region 1, both spinon and holon are in 
the 7r-flux state. The fluxes cancel completely. The region 3 
is the d-wave RVB phase where the d x 2_ y 2 pairing of spinon 
exists. The region 4 is the uniform RVB phase where hopping 
order parameters exist however it is real. In the region 5, all 
order parameters are zero. The phase diagram for t/J — 2 is 
quantitatively similar to the phase diagram for t/J = 1. 

We solved the self-consistency equations numerically, 
and obtained the phase diagram(Fig. 2) where 1) a re- 
gion of the staggcrcd-fhix phase of the electron exists, 2) 
staggered-flux and d-wave pairing do not coexist, and 3) 
the ground state is a purely superconducting state. 

At half-filling, the holon order parameter Bij is 
zero and the degeneracy of spinon states between the 
staggered-flux state and the d-wave pairing state ex- 
ists due to local SU(2) symmetry [34]; \ 7^ 0, B^ = 
0,2/g + rj 2 — const ^ 0. In region 1, the staggered-flux 
of electron does not exist although the staggered flux of 
spinon exists. The spinon and holon states are both ir- 
flux state [35] respectively. In the electron picture, the 
staggered-flux is canceled completely; \ 7^ 0, B ^ 0, <p s — 
4>h = tt,i] = 0. The d-density wave order parameter 
of the electron is y c = x-B sin (( 7r — 7r )/4) = 0. The 
staggered-flux order of the electron exists only in region 
2. The staggered current of electron exists in this phase. 
The spinon staggered-flux <j) s and holon staggered-flux 0h 
are not equal to n nor 0, and the holon staggered-flux am- 
plitude <j>h is not equal to the spinon staggered-flux am- 
plitude <j> 8 , <j) s ^ <j) h ; x ^ 0, B ^ 0, <j) s ^ 0, <p h ^ 0, i] = 0, 
and y e = xBs'm ((<j> s - 0h)/4) ^ 0. In region 3, d x i_ y i- 
wave pairing exists; \ ^ and B ^ 0, (f> s — cf>h = 0, r\ 7^ 0. 
and y c = 0. In region 4, there exists only uniform hop- 
ping order; \ ^ 0, B ^ 0, <j) s = (f> h = V = 0, and y e = 0. 
In region 5, all order parameters are zero. Spinon and 
holon cannot hop; x — B = cf> s = 0h = rj = 0. and 
y e = 0. The phase transitions, region 1 to region 3 and 



region 2 to region 3, are first order. Other phase transi- 
tions are second order. It is proved analytically that the 
staggered flux and the d-wave pairing do not coexist at 
finite doping. Details are given in Appendix. 

With the boson order parameter B^, the 7t-flux phase 
and staggered-flux phase of spinon and holon extends to 
higher-doped region compared to the previous work [23] , 
where B^ was not considered. The transition between 
region 1 and region 2 is a second order transition in our 
theory. (If one only focuses on spinon degree of free- 
dom, this does not look like phase transition [19,20,23].) 
There exists an order parameter that characterizes this 
transition. It is the staggered-flux of the electron. This 
is shown in Sec. IV B. 



Although the ground state is the d-wave superconduct- 
ing phase at finite doping, the staggered-flux phase exists 
at finite temperature. The staggered-flux state has larger 
entropy than the d-wave pairing state, due to the larger 
number of the excitation around finite size Fermi-surface. 
On the contrary, the Fermi-surface in the d-wave pair- 
ing state is always point-like. Therefore the staggered- 
flux state arises when temperature increases. These are 
shown in Sec. V. 



When the fluctuation around the saddle point solution 
(gauge field) is included, it is expected that the staggered 
flux phase extends to still higher doped region compared 
with the mean field phase diagram. There are three rea- 
sons, 1) the d-wave pairing of spinon will be destroyed, 

2) the instability to the staggered-flux state exists, and 

3) the symmetry breaking of the staggered-flux order is 
discrete. Let us elaborate these reasons. Firstly, the 
d-wave pairing of spinon is destroyed above the Bose- 
condensation temperature of holon when the fluctuation 
around the saddle point solution (gauge field) is included 
[36] . The d-wave RVB state without Bosc-condenstation 
does not occur, and the transition to the d-wave pairing 
state becomes the direct transition to the d-wave super- 
conducting state. The second reason is that the instabil- 
ity to the staggered-flux state has been discovered also in 
the gauge theory that includes the full lattice structure 
[37,38]. It is to the staggered-flux state at low doped re- 
gion and to the flux density-wave state (incommensurate 
staggered-flux state) at high doped region. The momen- 
tum that has anomaly is Q = (tt, it) to the staggered-flux 
state, and is (it, 7r — e) to the flux density-wave. At finite 
temperature, these anomalies occur. The third reason 
is that staggered flux state is more robust against the 
fluctuation than the d-wave pairing (d-wave RVB) state. 
This comes from the difference in the symmetry break- 
ings. The symmetry breaking in the staggered-flux state 
is discrete symmetry, Z-i. On the other hand, continuous 
symmetry, U(l), is broken in the d-wave pairing state. 
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B. Staggered-flux of the electron 

The doping dependences of the staggered-flux ampli- 
tudes are shown in Fig. 3. The electron staggered-flux 
amplitude (f> c is a difference of the staggered-flux ampli- 
tude of spinon (f> s and holon </>h Eq.(12). In Fig. 3 the 
phase changes as doping increases from region I to region 
2 and to finally region 3. In region I which is 7r-flux phase 
of spinon, staggered-fluxes of spinon and holon cancel 
completely and staggered-flux order of electron does not 
exist. However, in region 2, staggered flux amplitude is 
not 7r, fluxes does not cancel completely and staggered- 
flux amplitude of electron remains. The staggered-flux 
amplitude of holon decays rapidly as doping increases but 
it is not zero. At the transition point to d-wave pairing, 
staggred-flux order parameter becomes zero discontinu- 
ously. 



nl2 



Spinon flux 
Holon flux 
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S 

FIG. 3. The staggered flux order parameter of spinon <f> B , 
holon (f>h, and electron 4> c = 4>s — 4>h at T = 0.15J for t/J = 2. 
S is hole concentration. 



C. Spontaneous current of the electron 



Spinon current and holon current cancel locally, 



4 + 4 = o. 



(33) 



4 = i<(-g-Xy + tB^)flf ja - h.c), (34) 
Jt L j =it{x; j b\b j -h.c). (35) 

The (Xijfifj) term does not contribute to the current 

because it is real. Only (B*jfjfj) term contribute to the 
current. 

The local cancellation of spinon current and holon cur- 
rent does not exclude the possibility of the electron stag- 
gered current. In our formalism, the physical electron 
current is given by the imaginary part of the product of 
hopping order parameters of spinon Xij an d holon B*y 

4 = U ( c i c j ~ h.c.) 



= it{(flf j .)(b]b i )-c.c.} 

= 2t ImixijBtj} 

= 2t lm{xBe ±l( -^- M/i }. 



(36) 



The explicit form of the electron current in this sad- 
dle point is J?+ &i = (-iy2t X Bsm(cj) c /4:) = (-l) l 2ty c , 
J? +6i = -(-iy2txBsm(cj) c /A) = -(— l)*2tj/ e - The cur- 
rent of the electron is small at lower doped region because 
it is in proportion to the B = \(b\bj)\. At low doping, B 
is proportional to 5 and small. The magnetic field due to 
the orbital current of electron is same order with previous 
estimations of ICT 3 T [2,25]. 



The spinon current 4 an d the holon current 4 are given 

by 



V. FERMI SURFACE AND EXCITATION GAP 

AT (0,tt) 

The temperature dependences of Fermi surface and ex- 
citation gap at (0, 7r) are shown. These behaviors are con- 
sistent with angle-resolved photoemission spectroscopy 
(ARPES) experiments [7-14]. 



A. Fermi surface 



Fermi-surface is the loci of gapless excitation in mo- 
mentum space. It is defined by the zero-energy line 
of the electron excitation energy. The Fermi surface 
for 5 = 0.02 is shown in Fig. 4 at several tempera- 
tures. When temperature decreases, the Fermi-surface 
becomes smaller in the staggered-flux phase, and changes 
to a point-like Fermi-surface discontinuously at the tran- 
sition temperature to the d-wave pairing phase. How- 
ever, this discontinuity is small because Fermi surface in 
the staggered-flux phase has shrunk almost to a point at 
the transition point. The shape of Fermi-surfacc in the 
staggered-flux state of our theory is segment-like because 
the spectral weight has large enough value only along a 
part of the contour. 



G 



(a) 



t/J=2, 8=0.02 



ky/n 



0.5 



T=0.178J 
T=0.158J 
T=0.136J 
T=0 



\ 



0.5 

kx/K 



t/J=2, 8=0.02 




0.55 



FIG. 4. The temperature dependence 

of the zero-energy-line of fermion at S = 0.02. In (a), the 
first quadrant of the Brilluin zone is shown. In (b), expanded 
view around (7r/2,7r/2) is shown. In this doping, transition 
temperature to the d-wave pairing phase is about T = 0.135 J. 



The contribution of holon Green function comes mainly 
from the bottom of the holon band. The momenta at 
the bottom of the boson band are (0, 0) and (ir, ir) in 
the staggered flux state, and are (0, 0), (ir, 0), (0, n), and 
(ir, 7t) in the 7r-flux state. We can treat the bosons as if 
they are Bose condensed at these band minima. In this 
approximation, the spectral function of electron has the 
following form, 



A>, k) = -K% - E+) + M 2 ^ - ££)} 



where 



Hi 



(1 + e k /(4 + w k 2 )V 2 )/2, M 2 



(38) 
(1 



ekAe 2 . + W£) 1 / 2 )/2, A rn (uj,k) is the incoherent part of 



the spectral function, = ± y 'e 2 . + W£ 



trums of spinon in the staggered-flux state. In real- 
ity, the boson excitations exist at finite temperature. 
They make the discrete <5-function peaks have a finite 
width. This effect, however, will not change the shape 
of the Fermi-surface. Thus, we show the zero-energy- 
line of spi non in t he staggered-flux phase defined by 
E k = -\Ak + w k - M = in Fig. 4. The upper band 
spectrum E£ is always positive at finite doping where 
—fj, > 0. 

Although the zero-energy-line of fermion forms an el- 
lipse, the Fermi-surface can be considered as segment- 
like. The reason is that the intensity of spectrum func- 
tion ||uk| 2 is not symmetric between in the inner re- 
gion (where \k x \ + \k y \ < 7r) and in the outer region (where 
\k x \ + \k y \ > 7r). As shown in Fig. 5 it is stronger in the 
inner region. For example, on the line k x — k y , where 
d x 2_ y 2-ga,p due to staggered-flux order(W^) is always 
zero, the intensity is | at the inner side of the ellipse 
as |«k| 2 = f , but the intensity is zero at the outer side of 
the ellipse as |«k| 2 = 0. 



This result is quite different from the SU(2) mean-field 
theory, where the shape of Fermi-surface is always point- 
like because the particle hole symmetry(i.e.SU(2) sym- 
metry) of fermion remains even at finite doping. As the 
filling of fermions are always half-filled, the staggered 
flux state and d-wave pairing state can be assumed to 
be gauge equivalent. Artificial introduction of the phe- 
nomenological interactions are needed for the fermion to 
have a segment-like Fermi surface in the SU(2) theory. 

Now we give the details. The electron Green function 
G c is described by the product of spinon Green function 
G s and holon Green function G h [39]; 



t/J=2, 8=0.02 
T=0.178J 



0.45 




0.50 



krM 



0.55 



0.55 



, FIG. 5. The spectral weight A e (uj,\s) at uj = divided by 

G c {iLU ni k) = -TToY^ G S im,P)G h {ipe - «W„,p- k). (37) | (i.e. K| 2 ) at T = 0.178J. View from the (tt, 7r)-direction 
^ P „•„. „ is shown. 



Wt,P 



7 



B. Excitation gap at (0, n) 



VI. SUMMARY 



The temperature dependence of the excitation gap at 
(0, 7r) looks almost continuous at the transition point be- 
tween the staggered-flux phase and the d-wave pairing 
phase (Fig. 6), although it is a first order transition in 
the present mean field analysis. The energy at (0, 7r) are 
£sta gg crcd,±(0, 7r) = ± \4t yh + 3 Jy s /2 \ - ( i in the stagger ed 
flux state, and £d- P airin g ,±(0, ir) = ±\J \J? + (3Jr)/2) 2 in 
the d-wave RVB state, respectively. 



t/J=2 



-0.25J 
-0.5J 





1 1 

uniform " 




staggered flux 




-\ d-wave pairing - 




i = 



0.05 0.1 0.15 0.2 0.25 

T ( units of J ) 

FIG. 6. The temperature dependence of the excitation gap 
at (0, 7r) at 8 — 0.02 for t/J = 2. The phase changes as 
temperature decreases from the uniform phase to the stag- 
gered-flux phase and finally to the d-wave pairing phase. 



C. Comparison with experiments 



The ARPES experiments on Eh^S^CaC^Os+a 
(BSCCO) [7-13] and recent one on La2-5Sr,jCu04 
(LSCO) [14] showed that a segment-like Fermi-surface 
exists near the (ir/2,ir/2) in pseudo-gap region. [7-9]. 
There is the gap whose symmetry is d x 2_ y 2. The gap 
is largest near (ir, 0) and does not exists on the (0, 0)- 
(n, 7r) line. The reason that these Fermi-surface often 
called segment-like, or arc, is due to the shape where 
the spectral weight is finite. It does not form a closed 
curve. When one scan on the (0, 0)-(ir, 7r) line, the 
spectral weight only exists in inner region of Brillouin 
zone (where \k x \ + \ky | < tt) [8]. Thus, this is called 
segment-like. When temperature decrease, the segment- 
like Fermi-surface become smaller in pseudo-gap phase 
and change to the Fermi-point in the rf x 2_ y 2-wave super- 
conducting state [10,13]. On the other hand, the tem- 
perature dependence of the excitation gap at (0, tt) al- 
most looks continuous at the transition point between the 
pseudo-gap phase and the d x 2_ y 2--wave superconducting 
phase [13]. Our present results are consistent with these 
observations. 



The competition between the staggered flux, or the d- 
density wave, and the d-wave pairing is analyzed in two- 
dimensional t-J model based on U(l) slave boson mean- 
field theory. Not only staggered-flux of spinon but also 
staggered-flux of holon are considered, independently. In 
this formalism, the hopping order parameter of electron 
is described by a product of hopping order parameters 
of spinon and holon. The staggered-flux amplitude of 
the electron is a difference of staggered-flux amplitudes 
of spinon and holon. A phase diagram is obtained where 
1) a region of the electron staggered-flux state exists, 2) 
the staggered flux and the d-wave pairing do not coexist, 
and 3) the ground state is a purely d-wave superconduct- 
ing state. In 7r-flux phase of spinon, staggered-fluxes of 
spinon and holon cancel completely and the staggered- 
flux order of the electron does not exist. However, in 
the staggered-flux phase of spinon whose staggered-flux 
is not tt, fluxes do not cancel completely and staggered- 
flux order of electron exist. Thus, the phase transition 
between these two phases, 7r-flux phase and staggered- 
flux phase of spinon, becomes a second order transition 
in electron picture. The order parameter that character- 
izes this transition is the staggered-flux order parameter 
of the electron. The relation between the staggered cur- 
rent of electron and that of spinon is provided. The lo- 
cal cancellation of spinon current and holon current does 
not exclude the possibility of the electron staggered cur- 
rent. It is proved analytically that the staggered-flux and 
the d-wave pairing does not coexist except at half-filling. 
The condition for coexistence of staggered-flux and d- 
wave pairing as a saddle-point solution is provided. The 
instability of these two states, staggered flux and d-wave 
pairing, are discussed. The temperature dependence of 
following two quantities, Fermi surface and the excita- 
tion gap at (0,7r), are shown. The behaviors are consis- 
tent with ARPES experiments. When temperature de- 
creases, the segments-like Fermi-surface in the staggered 
flux phase becomes smaller and changes to a point at 
the transition temperature to the d-wave pairing phase. 
The temperature dependence of the excitation gap at 
(0, tt) looks continuous at the transition point between 
the staggered-flux phase and the d-wave pairing phase 
although these two phases are qualitatively different. 
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APPENDIX: COEXISTENCE CONDITION AND 
INSTABILITY 

It is proved analytically that the staggered flux and 
the d-wave pairing coexist only at half-filling. The in- 
stability of these two states, staggered flux and d-wave 
pairing, are also discussed. 

From self-consistent equations (27) (28) (30), following 
equation is derived. It must be satisfied for any states 
whose free energy is at saddle point. 







VsV 



(2i) 2 C(A+B) + ^B 



where the quantities, A, B, and C, are given by, 



-Y, 



7k- tanh(/?£ ks /2) 
2 -Eks 



tanh(/?£ ks /2) 



(s) 



c = -T- 

N 2(e k 2 + fff) 1 / 2 e P E L - 1 ' 



(Al) 

(A2) 
(A3) 
(A4) 



k,5 



where 7k- = cosk x — cosk y , A > 0, B < 0, C < 0. Here, 
the quantity B is different from B which describes the 
amplitude of the boson-hopping. B is zero at half-filling 
as /i = 0, and is negative at finite doping. C is zero at 
half- filling or when all holon condense at (0,0) or (7r,7r), 
and is negative when excited holon exists or when holon 
condense at (it, 0) or (0, ir) which is the situation of 7r-flux 
state. If the staggered flux and the d-wave pairing coex- 
ists, y s r] is non-zero. Then in order for the coexistence, 
eq.(Al) requires that the quantity in the square bracket 
must vanish. At half-filling, B = C = 0, so the condition 
is satisfied. At finite doping B < 0, so if C = there 
is no coexistence. When C ^ 0, we cannot expect the 
quantity in the square bracket vanish except at a special 
point in the phase space. However, from the continuity 
of the order parameters, it is impossible to have both y s 
and r\ nonzero only at a single point. 

The curvature of free energy at the saddle points of 
pure staggered-flux state and pure d-wave pairing state 
have the following forms respectively, 



[40]): In this case, pure d-wave pairing are always sta- 
ble, and pure staggered-flux state is unstable against 
infinitesimal d-wave pairing. The curvature at each 
state are (d 2 F ' / dy 2 )\ n ^. Va=0 = -2N(3J/4) 2 B > and 
(d 2 F/dT 1 2 )\ y ^o^o = 2A(3J/4) 2 B < 0. (iii) The case 
that B and C are both negative: A region where pure d- 
density wave is stable and a region where d-wave pairing 
is stable can both exist. 

The situation where B is zero and C is negative does 
not occur. B is zero only at half-filling where C is zero. 



g|^ a =o = 2A{-(^) 2 B + (2i) 2 C}, 
lv.*>,„=o = 2iV[(— ) B + — {l 



drj 2 



I 



4(2t)- 
3.7 



(A5) 
;}]■ (A6) 



The stability of the states depends on B and C. There 
are three cases: (i)The case that B and C are both zero: 
This is realized at half-filling, and (d 2 F '/ dy 2 )\ v ^o y =o = 
{d 2 F/dr] 2 )\ Vs ^o,r;=o = 0. We need not discuss this 
case. (ii) The case that B is negative and C is 
zero (this is the same situation that Zhang discussed 
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